Size-dependent enhancement of superconductivity in nanowires 
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A shape-dependent superconducting resonance can be expected when an energy level associated 
with the transverse motion in a wire passes through the Fermi surface. We show that the recently 
observed width-dependent increase of T c in Al and Sn nanowires is a consequence of this shape 
resonance effect. 
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Increasing the critical temperature (T c ) of a supercon- 
ductor (SC) has been a major challenge. On the one 
hand one can look for different materials which exhibit 
a higher T c . Such a search has been very successful over 
the last 20 years. On the other hand microstructuring 
of a superconductor is a different and new road which 
is able to modify T c (i.e. increase and/or decrease) and 
may also give us further insight in the basic mechanism 
of superconductivity. 

In earlier works on microstructuring of SC in the meso- 
scopic regime, enhancement of the critical current (j c ) 
was demonstrated to occur due to trapping of vortices. 
Also a large increase of the critical magnetic field (H c ) 
was realized through such mesoscopic structuring which 
is mostly a consequence of surface superconductivity. 
But in both cases the zero-magnetic field critical tem- 
perature was unaltered. The enhancement of j c and 
H c could be accurately described by phenomenological 
theories such as the London approach and the (time- 
dependent) Ginzburg-Landau theory. 

In the present Letter we are interested in modifying 
SC on the nano-scale. We deal with systems where 
the electron motion is limited to quasi-one-dimensions 
(ID). During the last decade nanowires have attracted 
much attention in the context of phase fluctuations of 
the order parameter (i.e. quantum phase slips) But 
quantization of the electron motion in the transverse di- 
rection was not investigated in much detail. However, 
very recently numerical investigation of the Bogoliubov- 
de Gennes equations has shown [2j that this quantiza- 
tion results in significant shape-dependent superconduct- 
ing resonances with a profound effect on the nanowire 
T c . Such systems have been the subject of recent exper- 
imental studies [1, 0, HI , and we demonstrate here that 
the width- dependent increase of T c found in these exper- 
iments is a manifestation of these shape resonances. 

More than 40 years ago, Blatt and Thompson Q cal- 
culated a remarkable sequence of peaks in the thick- 
ness dependence of the energy-gap parameter of single- 
crystalline superconducting nanofilms in the clean limit. 
They called these spikes shape resonances. At that time 
it was not possible to produce high-quality SCs with 
nanoscale dimensions (only very recently the thickness- 
dependent oscillations of T c were observed experimen- 



tally in ultrathin Pb films [7]). For decades atomic nuclei 
were the only systems where the interplay between quan- 
tum confinement and pairing of fermions could be studied 
experimentally and where the expectations of Blatt and 
Thompson were confirmed as a series of size resonances in 
the pairing energy gap of nuclei Q. Very recently high- 
quality nanowires have become available, where this res- 
onance effect is expected [|| to be significantly enhanced 
over the two-dimensional (2D) case. 

The physics of the shape-resonance effect in nanowires 
can be understood as follows. The superconducting or- 
der parameter is not simply the wave function of an or- 
dinary bound state of two fermions but the wave func- 
tion of a bound fermion pair in a medium. For example, 
in the homogeneous situation the Fourier components of 
the Cooper-pair wave function are suppressed for wave 
numbers less than the Fermi wave number due to the 
presence of the surrounding Fermi sea Q. Therefore, 
the Fourier transform of the Cooper-pair wave function 
is essentially nonzero only in the vicinity of the Fermi 
wave number. In general, for both homogeneous and in- 
homogeneous situations, the superconducting order pa- 
rameter is strongly dependent on No, the number of 
the single-electron states (for one spin projection) sit- 
uated in the Debye "window" around the Fermi level 
\p— fkOr), fi+fkur)], where [i is the chemical potential (th e 
Fermi energy) and is the Debye frequency [l(| fllj . 
More precisely, the mean energy density of these states 
per volume unit No/i^hiurjV) is of importance. When a 
SC is small enough, then this density varies together with 
its characteristic size. In particular, in the presence of 
quantum confinement the band of single-electron states 
in a clean single-crystalline wire is divided up in a series 
of parabolic one-dimensional subbands. While the width 
(the cross section a) of a wire increases/decreases, these 
subbands move down/up in energy. Each time when the 
bottom of a subband enters the Debye window, the den- 
sity N d / (2hjj dV) increases. This results in a sequence 
of peaks in the density of states versus a, i.e., the shape 
resonances. The resonances are significant for nanowires 
but are smoothed out with increasing nanowire width. In 
the large-cr limit the density No/i^huJrjV) is nearly con- 
stant and slowly approaching the well-known bulk value 
N(Q) = mkp I '(2iT 2 h 2 ) (see, for example, Ref. [ll|), where 
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fcp is the 31? Fermi wave number. It is of importance that 
the shape resonances in the density of states are accom- 
panied by significant peaks in the a— dependence of the 
averaged (over spatial coordinates) order parameter 
This average is close to the energy gap in the quasipar- 
ticle spectrum, whereas the latter (taken at zero tem- 
perature) is proportional to T c . Hence, the critical SC 
temperature in a nanowire is an oscillating function of 
a in the presence of shape resonances (these oscillations 
will be smoothened by fluctuations in the cross section). 

To explore shape resonances, one needs to deal with a 
single-crystalline SC (or with a polycrystalline SC made 
of strongly coupled grains) where the electron mean free 
path is comparable with the SC characteristic size. In the 
considered nanowires this path is at least of the order of 
the wire width [4, 5]. Such a nanowire is in the clean limit 
in the transverse direction: impurities influence only the 
electron motion in the longitudinal direction. Use of 
this approximation to geth er with arguments similar to 
the Anderson theorem [Kj makes it possible to conclude 
that nonmagnetic impurities do not produce a signifi- 
cant effect on the order parameter in this regime. The 
Bogoliubov-de Gennes (BdG) equations are the mi- 
croscopic equations describing a clean SC when the order 
parameter A(r) changes with position. In the absence of 
a magnetic field, A(r) can be chosen as a real quantity 
(phase effects are beyond the scope of the present work), 
and the BdG equations become 

Em(r) = (-^V 2 - ji)tti(r) + A(rH(r), (1) 

E lVl {v) = A(r)«,(r) - f-^V 2 - pWr), (2) 

where E$ is the quasiparticle spectrum, fi stands for 
the chemical potential and m denotes the electron band 
mass. Equations U) and ([2]) are supplemented by the 
self-consistency relation A = g J2i u i v i(^ ~ %fi) w 'rth. g 
the coupling constant and fi — f(Ei) the Fermi function. 
The summation is over all eigenstates with the kinetic en- 
ergy within the Debye window [p, — hujo, M + ?wd]- This 
kinetic energy is very close to the single-electron energy 
in a gas of unpaired electrons. The chemical potential is 
determined by n = (2/V) / d 3 r[\ Ui \ 2 h + h| 2 (l - /„•)] 
with n the mean electron density. Due to the quantum 
confinement in the transverse directions we have to set 
Ui{Y)\res — Vi(r)\ re s = on the wire surface while in 
the longitudinal direction periodic boundary conditions 
are used. 

Fig. [T] shows the numerical results of Eqs. (TT]) and © 
for a cylindrical wire with length L = 1/xm for Alparam- 
eters: gN(0) = 0.18 and hu D /k B = 375 K [liJGj. The 
bulk chemical potential is chosen as /itmik = 0.9 eV (it 
corresponds to n ~ 3.88 • 10 21 cm -3 , see the remark about 
influence of the electron-density choice below) and the 
electron band mass m is taken to be the free electron 
mass. In Fig. [T] (a) the zero-temperature energy gap A a 
normalized by its bulk value Abuik is plotted, whereas in 
panel (b) the mean density of the single-electron states 
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FIG. 1: (Color online) The width-dependent relative gap 
Acr/Abuik (a) and the relative mean density of the single- 
electron states in the Debye window Nd/C2^ujdVN(0)) (b) 
versus ct 1//2 for cylindrical Al nanowires with cross section a 
at zero temperature. 

in the Debye window No/(2hiUDV) is shown versus the 
effective nanowire diameter er 1 / 2 in units of N(0). Notice 
that shape resonances in No/ (2Hu>dVN(0)) are accom- 
panied by large oscillations in A CT / Abuik- Amplitudes 
of these oscillations are an order of magnitude larger as 
compared to similar peaks in the gap function of thin 
films with the same width @, [13]. We can point out, at 
least, three reasons for such enhancement. First, quan- 
tum confinement is much stronger in nanowires. Sec- 
ond, contrary to a film, most single-electron subbands of 
a nanowire are degenerate. For example, in the situa- 
tion of a cylindrical nanowire there are three quantum 
numbers for Ui and vc i = {j,m v ,k} with j the radial 
quantum number, m v the azimuthal quantum number, 
and k the wave number of the free longitudinal motion. 
Any energy subband with m v ^ is degenerate due to 
Ej jmv ,k = Ej,- mif ^- Third, even if the subbands are 
not degenerate, they can be situated very close to one 
another in energy (in particular, this is typical for sub- 
bands with large quantum numbers). In this case a res- 
onance has not decayed when the next one appears, and 
they join in one profound peak. Note that in thin films 
the superconducting shape resonances are well separated 
from each other and equidistant [1, [12] which is clearly 
not the case here. Another important point about Fig. [T] 
is that there is an important difference between depen- 
dencies presented in panels (a) and (b) of Fig. [TJ The 
energy gap is strongly enhanced at the resonant points 
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FIG. 2: (Color online) Critical temperature Tc.cr/Tc.buik ver- 
sus the square root of the cross section of a cylindrical Al 
nanowire. Symbols are the experimental results from Refs. 
0> HL ' the solid curve is a guide to the eye and indicates 
the general trend of the experimental results. The inset: the 
same dependence but in the small-cr region. 
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FIG. 3: (Color online) The shape resonances in the relative 
transition temperature for Al and Sn cylindrical nanowires. 



(with respect to the bulk value) while there are much 
less important drops between the resonances. This is dif- 
ferent from the density of states where the oscillations 
are almost centered around the bulk value. The reason 
is that the resonant increase of Nr>/ (2haj£)V) occurs due 
to states making practically the same contribution to the 
order parameter, and this collective mechanism plays the 
role of an " amplifier" (this is why the word " resonance" 
is appropriate here). Indeed, we have it.; oc e lkz and 
Vi oc e lkz , where z is the longitudinal coordinate. Hence, 
the product UiV* appearing in the self-consistency rela- 
tion for the order parameter does not depend on k and 
remains the same for a given subband. But this product 
strongly oscillates when the quantum numbers describ- 
ing the transverse states are changed which destroys the 
collective behavior beyond a resonant point. 

The resonant peaks in the energy-gap function are 
much higher than the bulk value, whereas drops of A CT 



between resonances are not so significant. The same is 
true for the critical temperature. This is the reason why 
the experimentally measured averaged critical tempera- 
ture increases monotonically with decreasing wire width 
in the presence of inevitable cross section variations in 
real samples. 

Fig. [2] shows the effect of the shape resonances on 
the transition temperature T C (T /T Cj buik for cylindrical Al 
nanowires. We calculated T c for discrete values of a 1 / 2 
with steps of O.Olnm. The nanowire length and the mean 
electron density are the same as in Fig. [TJ The the- 
oretical results are compared with recent experimental 
data for Al polycrystalline nanowires with strongly cou- 
pled grains. The three sets of experimental results corre- 
spond to three different initial samples having the trans- 
verse dimensions progressively reduced by the ion-beam 
sputtering method [jj. The squares @ and circles [l3| 
represent the observations made for two samples fabri- 
cated on the same chip, whereas the data given with the 
stars [ij are from another experiment. The horizontal 
bars are due to a distribution in the wire cross section. 
The vertical error bars are a consequence of uncertain- 
ties in the choice of the bulk superconducting tempera- 
ture due to the presence of impurities: the upper limit 
corresponds to the transition temperature in clean alu- 
minum X^buik = 1.19 K, the lower limit T c .buik = 1.26 K 
is used for the squares and circles (this is the transition 
temperature in the nanowire with a 1 ' 2 = 71 nm, the 
widest in these two sets) and T Cj b u ik = 1.27 K for the stars 
(T C (J for the nanowire with a 1 / 2 = 65 nm, the widest in 
the set of Ref. pj). The experimental transition tem- 
peratures shown in Fig. [2] were extracted from the tem- 
perature dependent resistance as the point at which the 
resistance dropped to one half of its normal-state value. 
From Fig. [2] it follows that the general trend of the exper- 
imental results is in good agreement with the calculated 
resonance amplitudes and is close to the average trend of 
resonances. 

The same conclusion holds for nearly single-crystalline 
Sn nanowires fabricated with the method of electrode- 
position into porous membranes [|[. Though there is no 
detailed information available about T c tT in this case, we 
have at our disposal a clear signature of its increase up to 
1.1 T c .buik in a cylindrical tin nanowire with a width of 
about 20 nm (tr 1 / 2 « 17.7 nm) [3]. Using the experimen- 
tal trend from Fig. [5] one obtains T Cj(T /r c ,buik ~ 1-25 for 
an aluminum nanowire with the same width. Hence, the 
size-dependent effect on the critical temperature is a fac- 
tor of 2 smaller when passing from Al to Sn for a nanowire 
with an effective diameter of about a 1 / 2 ~ 17.7nm. In 
Fig. [3] the transition temperature calculated with the 
BdG equations for a cylindrical nanowire is shown for Al 
and Sn (gN{0) = 0.25 and huj D /k B = 195 K 0, EH) 
near a 1 ! 2 — 17.7 nm. One can see that the ampli- 
tudes for Sn are indeed smaller, and the reduction factor 
changes between 1.4 and 2.5 (depending on the specific 
resonance) . 

These findings show that the shape-resonance effect 
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plays an essential role in the size-dependent increase of 
the critical SC temperature recently observed in experi- 
ments with Al and Sn nanowires. Note that the resonant 
amplitudes are weakly dependent on the total electron 
density in the metallic domain n « 10 21 — 10 23 cm~ 3 . For 
example, when the density rises from 3.88 • 10 21 cm~ 3 to 
2- 10 22 cm~ 3 , the amplitudes of the most profound shape 
resonances are reduced by 10% — 15% while the change in 
the "ordinary" resonant deviations from the bulk value 
is practically negligible 0]. The mean distance between 
two neighboring resonances is determined by the inverse 
3D Fermi wave number and is therefore very sensitive to 
the electron density n. For instance, instead of 1 — 2 res- 
onances per 0.2 nm for n — 3.88 • 10 21 cm~ 3 , we get 3 — 4 
resonances for n = 2 • 10 22 cm~ 3 Thus, the concrete 
value for n in the metallic domain is not very important 
for the main result of the present work. 



Concluding, we have shown that the size-dependent in- 
crease of the superconducting temperature recently found 
in Al and Sn nanowires is well explained by the shape- 
resonance effect. The surface-phonon mechanism |14j . 
which is usually considered to be responsible for the T c - 
enhancement in granular films and wires (note that this 
enhancement is insensitive to the cross section!) [l5l ]. can 
also contribute but in our systems it is of secondary im- 
portance. For example, for Al and Sn films with thickness 
of about 5 nm this contribution can be estimated to be 
rs 0.1K (l4| . While the amplitudes of the shape reso- 
nances are about IK for this film thickness (see Refs. Q 
and 03). 
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